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SUMI-iARY 


A theoretical solution is given for the critical stress of 
thin— walled cylinders loaded in torsion. The results are presented 
in terms of a few simple fonmilas and curves which are appllcahle 
to a wide range of cylinder dimensions from very short cylinders of 
leirge radius to long cylinders of small radius. Theoretical 
results are found to he in somewhat better agreement with experi- 
mental results than previous theoretical work for the same range 
of cylinder dimensions. 


INTRODUCTION 


For most practical purposes the solution to the problem of 
the buckling of cylinders in torsion was given by DonneD. in an 
important contribution to shell theory published in 1933 (reference 1) 
The present paper, which gives a solution to the same problem, 
has two main objectives: first, to present a theoretical solution 

of somewhat improved accuracy; second, to help complete a series 
of papers treating the buckling strength of curved sheet from a 
unified viewpoint based on a method of analysis ebsentially 
equivalent to that of Donnell but considerably simpler. (See, 
for example, references 2 and 3.) 

The method of solution in the present paper, is that developed 
In reference 3 * The steps in the theoretical computations of the 
critical stress are contained in the appendix. The results are 
given in the form of nondlmenslonal curves and simple approximate 
formulas which follow these curves closely in the usual range of 
cylinder dimensions. 


I 


MCA TN No . I 3 H 


SYMBOLS 


J,m,n 

P 

r 

t 

u 


w 

X 

y 

D 

E 

L 

Q 


®n-» 


Integers 

arbitrary constant 
radius of cylinder 
thickness of cylinder wall 

axial component of d? splacement; positive in x-direction 

circumferential component of displacement; positive in 
y-direction 

radial component of displacement; positive outward 
axial coordinate of cylinder 
circumferential coordinate of cylindex' 

Et^ 


J2(l - 


8 p 


flexural stiffness of plate per unit length 

Young's modulus 
length of cylinder 

mathematical operator defined in appendix 

curvature parameter j^l — or 

coefficients of deflection functions 

critical sheer-stress coefficient appearing In 

2 

formle T^j. = kg 2;^ 

(n2 + p2)2 + 


L^t 




n^(n2 + 02)2 




deflection functions defined in appendix 
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^ = J 


\ 

T 


cr 


half wave length of hucklea in circuaferential direction 
Poisson’s ratio 
critical shear stress 


V + 2 


a“ 


7 ”^ in-TBrse of defined hy v”**" = w 


EESliLTS AlCD DISCUSSION 


The critical shear stresses for cylinders are obtained from 
the equation 

n% 

■^cr = 


L^t 


The values of for cylinders with either simply supported or 

clamped edges are given in the form of logarithmic plots in 
figure 1. The ordinate in this figiu'e is ilae critical shear- 
stress coefficient kg The abscissa is a curvature parameter Z 
which is given directly by the theory and involves the dimensions 
of the cylinder end Poisson's ratio. 

For very short cylinders the value of the shear-stress coef- 
ficient approaches the values for flat plates, 5*3^ when the edges 
are simply supported and C .93 when the edges are clanged. As Z 
increases kg also increases and the curves which defined kg 
are given approximately by strai^^t lines. For siiq)ly supported 
cylinders, ' ' ■ 

kg = 0.85 

For cylinders' with clamped edges. 


k 


s 


0.93 
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The range of validity of these formulae is approximately 

100 < Z < 10 . 

t2 

For the case of long cylinders the curves of figure 1 split 
into a series of curves depending upon the radius— thickness ratio. 

These ctirves, vhlch correspond to buckling of the cylinder into 
two circumferential waves (n = 2), depart from the straight lines 

at approximately Z = 10^ or approximately ^ = 3^ • Because 

the critical shear stress of a long cylinder is almost 
Independent of end conditions, the curves for different values 
of r/t apply both to cylinders with simply supported edges and 
to cylinders with clamped edges. These curves are probably some- 
what inaccurate, however, because one of the requirements for the 
validity of the simplified equation of equilibrium used is 
that n2 » 1, a calculation for long cylinders made by Schwerin 
and reported in reference 1 by Donnell suggests that all values 
corresponding to the curves given in the present paper for n = 2 
are slightly high. 

In figure 2 the results of the present paper are compered with 
those given by Donnell (reference 1) and Leggett (reference 4) . 

The present solution agrees quite closely with that of Donnell . 
except in the transition region between the horizontal part and 
the sloping straight-line part of the curves. In this region the 
present results are appreciably lees than those of Donnell 
(maximum deviation about 1? percent) but are in close agreement with 
Leggett *B results, which are limited to low values of Z. 

In figure 3 the present solution and that of Donnell for the 
critical shear, stress of simply supported cylinders are compared 
on the basis of agreement with test results obtained by a number 
of investigators. (See references!, 5, 6, and 7.) The curves 
giving the present solution are appreciably closer to the test 
points . More than 8o percent of the test points are within 20 percent 
of the values corresponding to the theoretical curve for simply 
supported cylinders given in the present paper, and all points 
are within 35 percent of veilues corresponding to the curve. 

In figure 4 the present solution for critical shear-stress 
coefficients of long cylinders which buckle into two half waves 
is given more fully than in figure 1 and is compared with test 
results of references 1 and 8. 
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The computed values from which the theoreticai curves presented 
In this paper were drawn ere given in tables 1 and 2. 


CONCLUDING EEMAEKS 


A theoretical solution is given for the buckling stress of 
thln-walled cylinders loaded in torsion. The results are applicable 
to a wide range of cylinder dimensions from very short cylinders 
of large radius to very long cylinders of small radius. The 
t'iieoretical results are found to be in somewhat better agreement with 
experimental results than previous theoretical work for the same 
range of cylinder dimensions. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., March 20, 19^7 
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APPEroiZ 

THEORETICAL SOLIUTON 


The critical shear etrees at which buckling occurs in a 
cylindrical ehell may be obtained by solving the equation of 
equilibritim. 

Equation of equilibrium .- The equation of equilibrium for 
a slightly buckled cylindrical shell under shear is (reference 3) 


)V^w 


Et, 




— + 2T(;rt 




W 


ox^ 


5y 


= 0 


( 1 ) 


where x is the axial direction and y the circ^unferentlal 
direction. The following figure shows the coordinate system 
used in the analysis: 



j 


1 

t 




111 ' 
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Dividing through equation (1) hy D gives 


+ 12^0-4 ^ + 2k 


( 2 ) 


where the dimensionless parameters Z and kg are defined hy 


and 




jt2D 


The equation of equilibrium may be represented by 


Qw = 0 


(3) 


where Q is defined by 


Q = 






Method of solution .— The equation of equilibrium may be solved 

by using the Galerkin method as outlined in reference 9. In 
applying this method, equation (3) is solved by expressing w in 
terms of an arbitrary number of functions (Vq, Vi, . . . Vj, Wo, 
Wl, . . ., Wj) that need not satisfy the equation but do satisfy 
the boundary conditions on w; thus let 

j- L 

buiWin (4) 

m=0 m=d 


w 


8 
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The coefficient Sp and bm are then detenained by th© e<juatlons 



2XnL 


Uo to 
2 \ n L 

0 uo 


VjjQw dx dy = 0 




T(^Qw dr dy = 0 


( 5 ) 


where 

n = 0^ 1» J 


The solutions given In the present paper satisfy the following 
conditions at the ends of the cylinder: 


For cylinders of short and medium, length with simply supported 

j 

eages w. = = v = 0 and u is unrestrained. For cylinders of 

short and medium length wi th clamped edges w = ^ = u = 0 and v is 
unresLrained. For long cylinders w = 0. (See references 2 and 3 . ) 


Solution for Cylinders of Short and Medium Length 

Simply supported edges .— A deflection fxuictlon for simply 
supported edges my be taken as the infinite series 


05 ■ -Op' 

v = 8ln2K\ % Bill SS » cos S in' Bln B2 (6) 

X / ■ . L ■ , X L ! 

m=l ... . m;=l ' , 


where X is the half wave length of -the buckles In the cir- 
cumferential direction. Equation ( 6 ) is equivalent to equation ( 4 ) 
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• V 


% = 


sin sin 
\ L 


> 


W. 


n 


cos Bin 22S 
V T. 


(T) 


Sulsstitutlon of' expressions (6) and (7) Into equations (5) end 
Integi’atlon over the . limits Indicated give 


n 


<n2 +. p2>a + 


>t'»(n2 + pS)® I 


-®M\ V 


mn 


/ ‘'jn 

Jf 1)2 — u)2 

m=l 


r' 


>(8) 


\ I (n2 + p2)2 + 


12Z^n^ 


8lCaP \ 


mn 


+ “^ ; aa-- ~— 0 


«^(n2 + p2)2 I n /_ n2 - m2 

m=l 


vhere 



n == 1, 2, 3, . . . 


and m ± n is odd. Equations (8) have a solution if the following 
determinant vanishes: 


10 


n=l 

n=2 

n=3 

n=h 

n=5 

n =6 




2 

0 


a 


'•e 


0 0 

0 0 

0 0 

X) . 0 




% 

0 

0 

0 


° 


^5 

0 

0 

0 

0 


0 

0 


° 


% 

0 

0 

0 

0 

0 


0 ^ ••• 


n=l 

0 

.2 

3 

0 

4 

'15 

0 

n=2 

2 

3 

0 

"1 

0 

2l 

n=3 

0 

6 

5 

0 

12 

7 

0 

II 

4 

15 

0 

12 

7 

0 

20 

9 

n=5 

0 

12 

21 

0 

20 

9 

0 

n=6 

6 

35 

0 

2 

3 

0 

22 

11 

m 

• 

V 

m 

» 

9 

• 

* 

• 

m 

• 

• 

• 

• 

V 

• 

• 

. 
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‘35 

0 

2 

■3 

0 

.20 

11 

0 


^1 
0 

.2 

3 

0 

‘l5 

0 - 

6 
35 


2 

3' 


D3 

0 

6 

^ 5 

-f ° 

0 -f 

0 

2 
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^5 
0 

m 

21 

f ° 


Jl 

15 

0 


0 ^ 


10 

21 

0 


.2D 

9 

0 


• ° ° 

. 0 7^^^ 0 

0 

. 0 0 

0 0 0 

. 0 0 0 


0 

ks •> 


0 

0 

0 




20 

9 


0 

.22 

11 


0 ^4 


0 

0 

0 

0 


® ^5 


i. 

35 

0 

2 

3 


22 

11 

0 


0 

0 

0 

0 

0 


0 


° ^ 


where 


Mji = -^ 

^ 83 


(„2 . p 2)2 * 


)t’t(n2 + p2)^J 


By rearranging rows and columns, the infinite determinant can he factored 
into the product of two infinite suhdeterminants which are equivalent 
to each other. The critical stress may then he obtained from the 
following equation: 
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®1 

h2 



®5 

^6 

• « . 

ag 

^3 


^5 ag . . . 

n=l 


2 

3 

0 

k 

13 

0 

*6 

35 

... 0 

0 

0 

0 

0 0 ... 

n=2 

2 

I 


6 

'3 

0 

,51 

0 

... 0 

0 

0 

0 

0 0 . .. 

n=3 

0 

6 

^3 

12 

7 

0 

2 

7 

... 0 

0 

0 

0 

0 0 . . . 

-ri- 

ll 

4 

15 

0 

12 

‘7 

S''-* 

20 

‘ 9 

0 

... 0 

0 

0 

'f 

0 

0 0 ... 

n=5 

0 

10 

"21 

0 

.20. 

9 


30 

IT 

... 0 

0 

0 

0 

0 0 ... 

n=6 

• 

6 

35 

• 

0 

• 

2 

3 

• 

0 

• 

30 

11 

« 

• 

... 0 
• 

0 

• 

0 

• 

0 

• 

0 0 . . . 

• • 

• 

n=l 

m 

0 

« 

0 

• 

0 

• 

0 

« 

0 

• 

0 

• 

••• 

• 

2 

3 

• - 

0 

• 

• 

15 

• * 
• • 

0 -i. 

0 35 

n=2 

0 

0 

0 

0 

0 

0 

... 4 


0 

0 

^ 0 








3 

ks 2 

5 

12 

'T 

21 

n=3 

0 

0 

0 

0 

0 

0 

... 0 

6 

5 

1^3 

0 ... 

3 

n=4 

0 

0 

0 

0 

0 

0 

V • • * - 

15 

0 

-¥ 


— ^ 0 

9 ^ • 

n=5 

0 

0 

0 

0 ■ 

0 

0 

... 0 

IQ 

21 

0 

f ^ ’IT ••• 

n=6 

0 

0 

0 

■ 0 

0 

0 

6 

• • • " 

35 

0 

2 

’3 

0 

‘ll ^ ••• 

* 

• 

* 

• 

• 

• 

• 

• 

• 

• 

* 

• 

• 

• • 

# 

1 • 

m 

• 

• 

« 

• 

• 

# 

• 

• 

• 

• 

• • 

• • 


= 0 


( 10 ) 


The first approximation, oibtalned from the second- order determinant, 
is given hy 


= (I J m2 


( 11 ) 
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The Becond approxlniatlon^ obtained from the third-order determinant^ 
is given hy 


k 2 


M1M2M3 


. (afn. 


( 12 ) 


The third approximation, obtained from the fourth-order determinant, 
iB given by 


XI * h } - '^=1 " \ 3^“'* * (15)^3 * (f)'«3’\ 


+ = 0 


(13) 


Each of these equations shows that for a selected value of the 
curvature parameter Z the critical buckling stress of a cylinder 
depends on the wave length. Since a structure buckles at the lowest 
stress at which Instability can occur, kg is minimized with respect 
to the wave length by substituting values of 3 into the equation 
until the minimum value of kg can be obtained from of kg 

against p. This procedure is permissible when P ^ ^ * that is, 

when the cylinder buckles into, more than two circumferentiel waves. 

For the limiting case of a cylinder buckling into two waves, 
see the section of the present appendix entitled "Solution for 
a Long Cylinder" which follows. 

Figure 5 (a) shows the convergence of the determinant for cylinders 
with simply supported edges. 

ClflTnped edges .— A procedure similar to that used for cylinders 
with simply supported edges -may be followed for cylinders with 
damped edges. The deflection function used is the following 
series ; 


sin \ a;n 
mR) 


COB cos 


(m + 2)rtX 



+ cos ^ ^ bjj^ 

mR 


nwx (m + 2 )nx 

COB — cos = 

Xi 1 j 


( 14 ) 


IJN No. 1344 


13 


&ch tera tf thlp 50rl©e satleflce the coiuUtlcn on v at the edges. 
The functions Vij end Wh or© now deilned as follows; 


Vjj a Bln ^ 

A* 


W„ a COS 

^ X 


COS ESS _ cog 

L L 


cos ESS _ 008 Xs •*• 


(15) 


where 


n a 2j ... 


When the same operations as those carried out for the case 
or simply supported edges are perfcarmed, the following simultaneous 
equations result: 

For n = 0, 


+ M 2 ) - a2M2 + ^^s /_ bm 


m=l,3>5 L 


ei 2 ^ (m -f 2)^ 


' m' 


1^-4 (m + 2)2 — 4 


For n a 1, 


ai(Mi + M3) - S3M3 + kg ^ hm 

m=p,2,4 








Lm2 — 1 in? — 9 


(m ^ (m + 2)2 

(m + 2)2 _ 1 (m + 2)2 - 9 
For n = 2 , 3, 4 ... , 


= 0 


an(Mn + Mn+2) - aji_2Mn - a„+2Mn+2 + ^ 

m=0 



m2 

- (n + 2)^ 


^ fa +2)^ + fm + 2)2 

(m + 2)2 - n2 (m + 2)2 - (n + 2)2 


Ih 
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vhere min is odd. ' 
For n = 0 , 


bo(2Mo+ M2) 




00 

I 




m2 (m + 2)2 

^n;; ^7*2)2 


0 


For n = 1 , 
l>l(Ml + M3) 


CO 


— ^ 3^3 


V 
’ z_. , 

m= 0 , 2,4 






JTF 


m 2 — 1 m^ — 9 


(m + 2)2 
(m +2)2-1 


(m + 2)2 - 9 _ 


= 0 


For n = 2 , 3 , 4 , 


^ri(Mn ■*■ Mn+2) ” ^n-2Mn ” ^n+2Mn+2 4 q / am 

’ n=0 


m2 




m‘ 


i2 — n2 m2 - (n + 2) ' 


= rm + 2^2 ^ fm + 2)2 

(m + 2)2 - n2 (m + 2)2 - (n + 2)2 


(16) 


vhere min is odd and 


«n 


7 T 

8p 



p2)2 * 

rt^(n2 + 32)2 


The Infinite determinant formed hy these equations can he rearranged 
so as to factor into the product of two determinants which are 
equivalent to each other. The vanishing of one of these determinants 
leads to the following equation (lijnited for convenience to the 
sixth order) ; 


I li 


r 



TN No. 


15 



^0 


ag 

^3 

®4 


n=0 

^(2 Mo+M 2) 



_ 64 
” 105 

0 

32 
” 315 

n=l 

32 

15 

^(Mi+M 3) 

_ 352 
105 

-^3 

22 

35 

0 

n=2 


- 3^2 
105 

^(MgH-M^) 

1472 

315 

-k^ 

1376 

1155 

n=3 

6k 

~ 105 

“^3 

1412 

315 


4l6o 

693 


n=4 

0 

35 

- 

_ 4160 
693 


9440 

1287' 

n=5 

_ JL 
315 

0 

1176 

1155 


9440 

1287 



The first approximation, ohtained from the second-order 
determinant, is given hy 

(2 Mq + Mg) (Ml + M3) (18) 


The second approximation, obtained from the third— order determinant, 
is given by 


(Ml + M3) [ (aip + Mg) (M2 + M4) - M2^] 

(ST il'a *{Mf 


(19) 


The third approximation, obtained from the fourth-order determinant, 
is given by . 
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cylinders with almply supported ed^es. Figure 5(b) shows the convergence ; of this 
determinant. 
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Solution for e Long Cylinder 

A long slender cylinder ^ 10 ^ will ‘buckle into two 

waves in the circumferentled direction. If, in the previous 
cases of cyllgders with 8in5)ly supported or clanped edges, the 
half wave length in the circumferential direction X is taken 
as «r/2, it is possible to find the critical stress of a long 
Blender cylinder with the corresponding edge conditions. This 
method of solution is laborious, however, because determinants of 
high order must be employed to obtain solutions of reasonable 
accuracy. The labor is greatly reduced by the use of the following 
deflection function: 


V. = 




(p + 2 )kx 
L 


+ 



(21) 


where p + 1 is the phase difference of the circumferential waves 
at the two ends of the cylinder meeisured in quai'teiwevolutlons. 
This equation satisfies the single boundary condition w = 0. 

With this deflection function, the functions V and V all 
vanish except 


Vl 


cos 


( 


+ 2y V (p_HLg),2g. + gy 

r / L ^ ^ 


(gg) 


Use of equations (5), (gl)^ and (22) and the relation 2X = rtr 
results in the following equation: 


r 



(P + g)? + 




12Z^(p H- 2) 


[_(P + + 


y ( 23 ) 


33 
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This equation may te written 



( 24 ) 


• For given values of Z and ^ j/l — P Is varied until a 

minimum value of ICg is obtained from a plot of p and corresponding 
values of kg. The critical stress of e long slender cylinder is 
very insensitive to edge restraint; therefore, the solution applies 
with sufficient accuracy to cylinders with either simply supported 
or clamped edges. The Ehear-stress coefficient for long slender 
cylinders is plotted against the curvature parameter in figure 4, 
and parts of these curves also appear in figure 1. 


I I 
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i TABLE 1 


THEORETICAL SHEAR-STRESS COEFFICIENTS AND WAVE LENGTHS 


OF BUCn^S FOR SHORT- AND MEDIUFt-LENGTH CYLINDERS 


i 

z 

i 

1 

First approximation 

Second approximation | 

i 

Third app: 

rozimation 


3 


3 

J^s 

P 

Cylinders with simply supported ed^ss 

0 

5.60 

0.770 

5.34 

0.790 



1 

5.69 

.805 

5.42 

. 36 o 

5.41 

0.865 

5 

6.68 

1.00 

6.22 

1.015 



10 

8.36 

1.24 

7.55 

1.265 

7.545 

1.27 

30 

14.93 

1.82 

12.69 

1.875 



100 

34.09 

2.74 

27.86 

2.91 



300 

76.80 

3.86 

62.47 

4.18 

61.47 

4.32 

1,000 

189.5 

5.4o 

153.0 

5.95 


10,000 

107? 

10.0 

871.2 

11.2 

851.9 1 

11.8 


6050 

17.9 

4920 

20.1 


23.0 

Cylinders with cla 

mp«d edges 

0 

9.55 

1.175 

9.31 

1.205 

9.09 

1.205 

1 

9.57 

1.18 

9-32 

1.21 



5 

9.90 

1.23 

9.62 

1.27 



10 

10.79 

1.35 

10.42 

1.38 

10.19 

1.38 

30 

16.13 

1.89 

14.99 

1.97 


• 100 

35.40 

2.95 

30.68 

3.14 

30.65 

3.12 

1,000 

206.3 

6.12 

167.5 

6.70 

165.7 

7.00 

10,000 

6860 

20.35 

5449 ' 

23.2 

5310 

24.8 
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Figure 1.- Critical shear-stress coefficients for thin-walled cylinders 

in torsion. , 
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